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We expand on our earlier work (cond-mat/0612130, Phys. Rev. Lett. 99, 46404 (2007) ) in 
which we constructed the exact low-energy theory of a doped Mott insulator by explicitly integrating 
(rather than projecting) out the degrees of freedom far away from the chemical potential. The exact 
low-energy theory contains degrees of freedom that cannot be obtained from projective schemes. In 
particular a new charge ±2e bosonic field emerges at low energies that is not made out of elemental 
excitations. Such a field accounts for dynamical spectral weight transfer across the Mott gap. 
At half-filling, we show that two such excitations emerge which play a crucial role in preserving 
the Luttinger surface along which the single-particle Green function vanishes. In addition, the 
interactions with the bosonic fields defeat the artificial local SU(2) symmetry that is present in 
the Heisenberg model. We also apply this method to the Anderson-U impurity and show that in 
addition to the Kondo interaction, bosonic degrees of freedom appear as well. Finally, we show that 
as a result of the bosonic degree of freedom, the electron at low energies is in a linear superposition 
of two excitations-one arising from the standard projection into the low-energy sector and the other 
from the binding of a hole and the boson. 



I. INTRODUCTION 

Low energy theories based on an explicit integra- 
tion over the degrees of freedom at high energy are 
the cornerstonei. for analyzing long wavelength physics 
of interacting systems. For high-temperature super- 
conductivity in the cuprates, the relevant^ low-energy 
theory must be constructed for a doped Mott in- 
sulator. While no shortage of theories has been 
proposed^iii^i^ili^'iaiiiii^^^^i^iili^iiaiSOiSi^ none is 

based on an explicit integration over the degrees of free- 
dom at high energy. The primary difficulties in carrying 
out such a program appear in the simplest model 

-ffHubb = 5y4aCi,<T -I- j^C,,J,Ci,| 

- Ht + Hu (1) 

applicable to a doped Mott insulator. Here i,j label lat- 
tice sites, gij is equal to one iff i,j are nearest neigh- 
bours, Ci,a- annihilates an electron with spin cr on lattice 
site z, t is the nearest-neighbour hopping matrix element 
and U the energy cost when two electrons doubly oc- 
cupy the same site. The cuprates live in the strongly 
coupled regime in which the interactions dominate as 
t w 0.5eV and U w 4eV. As U is the largest energy 
scale, it is appropriate to integrate over the fields that 
generate the physics on the U scale. The operators that 
correspond to such physics can be written in terms of 

??i,cr = Ci^aTli-a, noting that ni^^Ui^i = J2aVl,crVi,cr/2. 

Physically, rji^cr ivl a) annihilates (creates) an electron on 
a doubly (singly) occupied site, hence is associated with 
the energy scale U. Consequently, the interaction term 
reduces to a simple quadratic form 

Hu ^U^n,^^n,^l = ^^r]l^r],^„, (2) 



which would enable an exact integration over the U scale 
*/ Vi-u obeyed canonical fermionic or bosonic commuta- 
tion relations. However, a simple computation gives 

Hence, standard bosonic or fermionic coherent state rep- 
resentations are of no use in integrating over the fields 

The additional problem is spectral weight transfer. 
When one electron resides on each site (half-filling), a 
charge gap of order U opens for all t/ in d = 1 and 
provided U ^ t for d > The band above the 

gap describes electron motion on singly occupied sites, 
whereas the band below captures electron motion on 
empty sites. Such motion is described by rji^^r and 
^i,a — Cj.o-(l — ni,^a-), respectively. However, unlike the 
traditional band picture in which electron motion occurs 
in either the conduction or valence bands, electron spec- 
tral weight lives both above and below the Mott gap. 
This state of affairs obtains because the electron annihi- 
lation operator 

^i.cr — ^i,a{^ ^i, — cr) Ci,a^i, — a — ^i,a ^i.cr (4) 

can be written as a linear combination of excitations that 
reside in both bands. That is, unlike the standard band 
insulator picture, the states lying above the gap are not 
orthogonal to those below it. As a consequence, adding 
or removing electrons from a Mott insulator changes the 
distribution of spectral weight at all energies. In particu- 
lar, the addition of x holes to a Mott insulator creates at 
least 2x single particle addition states^'^ just above the 
chemical potential. The deviation from would be 

the case in a band insulator, is intrinsic to the strong 
correlations that mediate the Mott insulating state in 
a half-filled band, thereby distinguishing Mottness from 
ordering. Each hole reduces the number of ways of cre- 
ating a doubly occupied site by one, thereby reducing 
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the spectral weight at high energy. As the x empty sites 
can be occupied by either spin up or spin down electrons, 
the 2x sum rule is exaclj^^ in the atomic limit, U ^ oo. 
In the presence of hybridization (with matrix clement t) , 
virtual excitations between the LHB and UHB increase 
the loss of spectral weight at high energy thereby lead- 
ing to a faster than 2x growth^^-^^'^^ of the low-energy 
spectral weight, a phenomenon confirme d^^i^^'^^ widely 
in the high-temperature copper-oxide superconductors. 

A true low-energy is not exact if it cannot account for 
all low-energy degrees of freedom even if they arise from 
the high-energy scale. Hence, the true low-energy theory 
of a doped Mott insulator must preserve the 2x sum rule. 
In this regard, two approaches are possible: CI) change 
the particle statistics so that placing a particle on one site 
excludes particles of opposite spin or C2) generate new 
degrees of freedom at low energy so that removal of an 
electron destroys at least charge e states. Perturbative 
methods followed by projectio n^^'^'^i^^'^^i^^ of the high- 
energy scale as well as slav o^^'^°i^^ particle techniques all 
implement CI. To leading order in t^/J7, the result is the 
t— J model. The key goal in such approaches is to diago- 
nalize the Hubbard model into sectors with a fixed num- 
ber of doubly occupied sites. When one performs such a 
transformation^^, however, the electron operators must 
be transformed as well. Although this step is generally 

^„j^Qj.^,j3. 4. 5. 6. 7. 8. 9 . 10. 11. 12. 13. 14. 15. 16. 17. 18. 19. 20. 21 jg 

cial because the no double occupancy condition applies 
only to the transformed fermions not to the bare elec- 
trons. As the relationship between the transformed and 
bare electrons is non-linear, it is advantageous to devise 
a much simpler method in which the mixing to the dou- 
bly occupied sectors in the bare electrons is carried by a 
single degree of freedom. The current method provides a 
solution to this problem. All the physics associated with 
the mixing between the UV and IR scales is captured 
by a charge 2e bosonic field. In addition, one might en- 
tertain the possibility that slaved-particle method o^^'^°'^^ 
could be tailored to implement an integration of the high- 
energy scale. However, in the slaved operator approach, 
the interactions involving double occupancy are highly 
non-linear as a result of the constraints that remove the 
unphysical states and hence double occupancy cannot be 
integrated over explicitly. 

We present here a detailed description of a method that 
permits an explicit integration of the degrees of freedom 
far away from the chemical potential in doped Mott in- 
sulators. We show that the degrees of freedom far away 
from the chemical potential can be integrated out ex- 
plicitly (that is, without resorting to projection or slave- 
particles) for a doped Mott insulator. The result is that 
new physics emerges at low energies, namely a charge 2e 
boson, that cannot be thought of as simply related to 
electronic motion. Our work lays plain that the true-low 
energy theory of the Hubbard model is not at — J-like 
model in terms of the original bare electrons as is com- 
monly believed. The true low-energy theory is an exam- 
ple of C2. 



The charge 2e boson enters the theory as a Lagrange 
multiplier field. As such, it does not have dynamics, 
in much the same way that the a field in non-linear a- 
models does not have dynamics. In that case, dynamics 
are generated radiatively, by taking into account interac- 
tions with the other fields in the model. In the simplest 
spherical case, the latter fields can be completely inte- 
grated out, and a large TV expansion organizes the theory. 
The situation for Mottness is considerably more compli- 
cated. In particular, we have not yet elucidated the pre- 
cise low energy dynamics. Instead, we take the appear- 
ance of the 2e boson as an indication that the building 
blocks for the low energy dynamics of strongly correlated 
electron matter involve degrees of freedom that lack elec- 
tron quantum numbers. As we show in the present paper, 
and more fully in a companion papers22i^, there are in- 
dications that the boson should not be thought of as a 
weakly interacting dynamical field at low energies with, 
for example, a Fock space of its own, but that instead 
it should be thought of as a constituent in a strongly 
coupled theory. For example, from the exact form of 
the electron creation operator at low energy, we deduce 
that the boson can mediate new charge e excitations by 
binding a hole. It is the emergence of this state at low 
energies that serves to preserve the 2x sum rulc^"^. We 
believe that there are analogies here between the pres- 
ence of such composite states and confining dynamics in 
particle physics. Indeed, the nature of an insulator is of 
course that electric transport is absent, in analogy to the 
absence of color transport in QCD. 

This work expands considerably our previous Letter- 
in which we presented only an outline of the method. 



II. LOW-ENERGY THEORY 

We will be concerned with the limit in which the Hub- 
bard bands are well-separated, U ^ t. Given that the 
chemical potential lies in the gap between such well sep- 
arated bands at half-filling, which band we should asso- 
ciate with high energy is ambiguous at half-filling. Both 
double occupancy (UHB) and double holes (LHB) are 
equally costly. Doping removes this ambiguity. Hole- 
doping jumps the chemical potential to the top of the 
LHB thereby defining double occupancy to be the high 
energy scale. For electron doping, the chemical poten- 
tial lies at the bottom of the upper Hubbard band and it 
is the physics associated with double holes in the lower 
Hubbard band that must be coarse-grained. At half- 
filling, both the UHB and LHB must be integrated out. 
As each of these limits results in a different theory, we 
will present each separately. As will be seen, the low- 
energy theories that result from the electron and hole- 
doped cases are related, though not by the naive particle- 
hole transformation. 
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FIG. 1: (a) extended Hilbert space for a single site, (b) hop- 
ping processes between neighbouring sites included in the 
Lagrangian. Double occupation has been replaced by D- 
occupation. 



A. Hole Doping 

Within the Hilbert space for the Hubbard model, 
®i (J^l (8) .Fj), it is impossible to integrate out the de- 
grees of freedom far away from the chemical potential. 
The basic idea of our construction is to rewrite the Hub- 
bard model in such a way as to isolate the high energy 
degrees of freedom so that they can be simply integrated 
out. To solve this problem, it is expedient to extend 
the Hilbert space (^i [T^ ® Ti® Tb)- The key idea is to 
associate with the creation of double-occupation, to 
be implemented by a constraint. In order to limit the 
Hilbert space to single occupation in the D sector, we 
will take D to be fermionic. We refer to D as a fermionic 
oscillator as it is associated to a two state system. The 
field D will enter the theory as an elemental field with a 
large (order J7) quadratic term and precise interactions 
with the electronic degrees of freedom; the low-energy 
(IR) theory is obtained by integrating out D. The inter- 
actions of this extended model must be chosen so that the 
model is precisely equivalent to the Hubbard model; in- 
deed, if instead of integrating out the field D, we merely 
solve the aforementioned constraint, the model will re- 
duce to the Hubbard model, which we will refer to as the 
high energy (UV) theory. 

The action of the standard electron creation operator, 
and the new fermionic operator, D^, to create the 
allowed states on a single site are shown in Fig. [T] 

There are of course several unphysical states in this 
Hilbert space. As we will see, such states are removed 
once the constraint is solved. At present, the expan- 
sion of the Hilbert space should be thought of as a tool 
to enable the integration of the high energy degrees of 
freedom. To proceed, we formulate a Lagrangian in the 
extended Hilbert space. The allowed hops involving the 
D fields and the electron operators which are equiva- 
lent to the hops in the Hubbard model are indicated in 
FigHb). For example, the hopping process in the upper 
left-hand corner describes the hopping of a hole in the 



lower-Hubbard band. The terms in the middle describe 
the transport between the D field and two electrons in a 
singlet on neighbouring sites. The term in the lower right 
corner describes a hop in which Di and switch places. 
There are no further allowable hopping processes. A fur- 
ther requirement of the Lagrangian for the hole-doped 
theory is that it contain the appropriate dynamical term 
for motion in the lower Hubbard band. That is, those 
sites which contain the occupancy |c| ^ |0) must be ex- 
cluded from hopping processes (such hops are accounted 
for by the hopping of D). The Euclidean Lagrangian in 
the extended Hilbert space which describes the hopping 
processes detailed above can be written 



-I- (pyBci^^V^Cj.-^ + h.c.) 



(5) 



Here, gij selects out nearest neighbours (note that if we 
wanted to include next-to- nearest neighbour interactions, 
we need just modify the matrix gij accordingly), the pa- 
rameter Va has values V| = 1, V| = — 1, and simply en- 
sures that D couples to the spin singlet and the operator 



Cij,a is of the form C^. 



eoa, 



0(1 



n-j.-a) with number operators Ui^a- = cl^Ci^a-- 

For simplicity, we have introduced a complex Grass- 
man constant 9, which we have inserted in order to keep 
track of statistics; it bears some resemblance to a super- 
space coordinate. Because Dj is fermionic and Cj^i^Cj^i 
transforms as a boson, a Grassman variable is needed 
to essentially 'fermionize' double occupancy. They are 
normalized via 



= 1. 



(6) 



The Grassmann variable is an artificial device that will 
disappear from the UV or IR Lagrangians. 



The constraint Hamiltonian Hnon is taken to be 



ec 



h.c. 



(7) 



where is a complex charge 2e bosonic field which enters 
the theory as a Lagrange multiplier. The constant s has 
been inserted to carry the units of energy. At this point, 
there is some ambiguity in the normalization of ip, but 
we expect that this will be set dynamically. We will find 
that if a true infrared limit exists, then s must be of 
order the hopping matrix element t. There is a natural 
parallel between the constraint condition, Eq. ([7]), and 
the constraint in the non-linear sigma model. In fact, 
the auxiliary field ip will enter the low-energy theory in 
an analogous fashion to a in the non-linear sigma model. 
In both cases, if and a enter as Lagrange multipliers. 
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Both end up playing a crucial role in the phase structure 
of the true low-energy or infrared limit. In this case, ip 
will serve to create new excitations at low energy which 
will generate the dynamical part of the spectral weight 
transfer across the Mott gap. 

Now, as remarked previously, we have chosen the La- 
grangian (O so that this theory is equivalent to the Hub- 
bard model. To demonstrate this, we first show that once 
the constraint is solved, we obtain the Hubbard model. 
Hence, the Lagrangian we have formulated is the Hub- 
bard model written in a non-traditional form ~ in some 
sense, we have inserted unity into the Hubbard model 
path integral in a rather complicated fashion. To this 
end, we compute the partition function 

Z = j[Dc Vc^ VD VD^ Vip Vip^] exp" ^dt ^ (g) 

with L given by ([5]). We note that is a Lagrange mul- 
tiplier. As shown in the Appendix (Eq. (|4T|) ). in the Eu- 
clidean signature, the fluctuations of the real and imagi- 
nary parts of ifi must be integrated along the imaginary 
axis for tpi to be regarded as a Lagrangian multiplier. 
The if integrations (over the real and imaginary parts) 
are precisely a representation of (a series of) (5-functions 
of the form. 



dOD, 



(9) 



If we wish to recover the Hubbard model, we need only 
to integrate over 1?^, which is straightforward because of 
the (5-functions. The dynamical terms yield 



d^e ee 



^9 OO'Y^ [(1 - _<^)cJ_^Ci,o- -I- ni^-acla^i,o- 

J d^e eeY^i.ir,^- 



(10) 



Likewise the term proportional to yields 

hi 

= J d^9 86 ^ gy-rij _o.c]^Ci^cr -I- h.c. 



+ h.c. 



(11) 



Finally, the hopping terms that involve two D fields give 
rise to 

j d^eeoY^g^j [444|(4|C,,| -|-c]_^c,,4)cj-|Cj-4 



(12) 



Eqs. (fTTj) and (fT2|) add to the constrained hopping term 
in the Lagrangian (the term proportional to Cij^a) to 
yield the standard kinetic energy term in the Hubbard 
model. Finally, the D^D term generates the on-site re- 
pulsion of the Hubbard model. Consequently, by inte- 
grating over ipi followed by an integration over Di, we 
recover the Lagrangian, 



J d^e ^^LHubb = ' 



Hubb, 



(13) 



of the Hubbard model. This constitutes the ultra-violet 
(UV) limit of our theory. In this limit, it is clear that 
the Grassman variables amount to an insertion of unity 
and hence play no role. Further, in this limit the ex- 
tended Hilbert space contracts, unphysical states such as 
|1,0,1), |0,1,1), |1,1,1) are set to zero, and we identify 
|1,1,0) with |0,0, 1). Note there is no contradiction be- 
tween treating D as fermionic and the constraint in Eq. 
([7]) . The constraint never governs the commutation rela- 
tion for D. The value of D is determined by Eq. ([7]) only 
when (p is integrated over. This is followed immediately 
by an integration over D at which point D is eliminated 
from the theory. 

The advantage of our starting Lagrangian over the tra- 
ditional writing of the Hubbard model is that we are able 
to coarse grain the system cleanly for U ^ t. The en- 
ergy scale associated with D is the large on-site energy 
U. Hence, it makes sense, instead of solving the con- 
straint, to integrate out D. The resultant theory will 
contain explicitly the bosonic field, ip. As a result of this 
field, double occupancy will remain, though the energy 
cost will be shifted from the UV to the infrared (IR). 
Because the theory is Gaussian, the integration over Di 
can be done exactly. This is the ultimate utility of the 
expansion of the Hilbert space - we have isolated the 
high energy physics into this Gaussian field. As a result 
of the dynamical term in the action, integration over D 
will yield a theory that is frequency dependent. The fre- 
quency will enter in the combination uj + U which will 
appear in denominators. Since U is the largest energy 
scale, we expand in powers of w/U; the leading term 
yields the proper w = low-energy theory. Since the 
theory is Gaussian, it suffices to complete the square in 
the I?-field. To accomplish this, we define the matrix 



E 



(14) 



and bi = J2j ^ij = Y^j.a 9ij<^j.,a-VaCi-c. At zero frequency 
the Hamiltonian is 



5 



where 

i,j 

3 id 

which constitutes the true (IR) hniit as long as the energy 
scale s is not of order U. If s ^ 0{U) then we should 
also integrate out ipi - this integration is again Gaus- 
sian and can be done exactly; one can easily check that 
this leads precisely back to the UV theory, the Hubbard 
model.-^ Hence, the only way in which a low-energy the- 
ory of the Hubbard model exists is if the energy scale for 
the dynamics that ip mediates is 0{t). This observation 
is significant because it lays plain the principal condition 
for the existence of an IR limit of the Hubbard model. 

To fix s, we note that the fourth term entering our 
Hamiltonian can mediate spin exchange. As the energy 
scale for this process is <^/J7, we make the identification 
s ~ i. Hence, the low-energy theory contains a charge 
2e bosonic field which can either annihilate/create dou- 
bly occupied sites or nearest-neighbour singlets. That 
the energy cost for double occupancy in the IR is t^/J7 
and not U underscores the fact that the UHB and LHB 
are not orthogonal. The presence of the new field ipi at 
low energies is the result of the overlap between the high 
and low energy scales. Physically, double occupancy oc- 
curs at low energies for two distinct reasons. The first 
is spin exchange which is generated by the term (pibl- 
The second is motion of a doubly occupied site (a dou- 
blon) entirely in the LHB. The latter is not present in 
projective models but is certainly a low-energy process 
that must be present in the exact low-energy theory. 

While electron number conservation is broken in the 
IR, we find (by inspection of p3|) ) that a conserved low- 
energy charge does exist^i 

As Eq. (jlSp makes clear, bosons acquire dynamics only 
through electron motion. Hence, the low-energy theory 
of a hole-doped Mott insulator is a strongly coupled bose- 
fermi model. On purely phenomenological grounds, bose- 
fermi models have been advanced^i^ as a starting point 
for tackling the cuprate problem. In such models, and 
others^, the bosons are viewed as non-interacting and 
possess a Fock space of their own. However, the current 
analysis lays plain that while the bosonic degree of free- 
dom exists, it does not extend the Hilbert space of the 
Hubbard model. That is, the charge 2e boson does not 
have a Fock space of its own. In obtaining the low-energy 
theory, we integrated over the high energy D— field which 
acted in the extended Hilbert space. Consequently, the 
resultant low-energy theory preserves the original Hilbert 
space of the Hubbard model. As we will show, a distinct 
possibility is that the boson acts to create composite ex- 
citations that have charge e. 



Several limits are of interest. First, consider the limit 
U = oo (for fixed lattice size). The theory reduces to the 
restricted hopping term and the third term in Eq. (fT5|) . 
In this limit, the tp integration reduces to a delta function, 
i5(ci,|Ci^j^), giving a constraint enforcing the vanishing of 
double occupancy, the correct result for U = oo. 

Second, should (p = 0, we recover the interactions in 
the t—J model. To establish this, we note that for (pi = 0, 
we have the restricted hopping term and the first term in 
Eq. p^ . Approximating J\4ij by its leading term, Sij, 
the second term reduces to 

Y^l^^= 9i3 9Ucl^„Vacl^Ct,a'Va'C,-a' (17) 

i ijtaa' 

which contains the spin-spin interaction — (S.^ • Sj — 
niUj/A) as well as the three-site hopping term. Next, 
we expand the Tr In term 

TrhiM^^ = Tr\n{6ij + Aij) 

— An -f —AijAji + . . . 

where 

= Y 5y S>c»,o- (18) 

(7 

is nonzero only if i,j are nearest neighbours. When this 
term appears in the Euclidean Lagrangian, its magni- 
tude is kBTt^ /{to + Uy. Therefore, at low temperature, 
this term is small compared to and to leading or- 
der in t'^/U, the terms in Hi^ dominate. Hence, the 
ip — limit contains the interactions in the t—J model, 
thereby establishing that the physics contained in tpi is 
non-projective. To make closer contact with the t — J 
model in which the spin-spin interaction acts only in the 
singly-occupied sector, we note that the theory we have 
developed here could have been formulated strictly in the 
projected space by simply substituting ^i^a- for Ci^^ in the 
hopping terms containing Di in our starting Lagrangian, 
Eq. ([5]). The only substantive difference would be that 
the second hopping term (the term quadratic in Di) in 
Eq. ^ would enter with the opposite sign. Hence, in 
the IR hmit, Mij 6ij +t/U^j^^£,i^^. This change is dic- 
tated by the commutation relations of the ^^^o- operators. 
The UV limit, the Hubbard model, is obtained as before. 
Setting ip>i = Q va the IR limit leads exactly to the t—J 
model. Thus, the t — J modeliS written in terms of the 
bare electron operators is not the low-energy limit of the 
Hubbard model. This is not entirely surprising as Es- 
kes and others^^ has stressed that the operators must be 
transformed as well in writing the t — J model. Only at 
U = oo do the transformed and bare fermion operators 
agree. Hence, at any finite U , the physics is governed 
by a finite length scale for double occupancy. Hence the 
limits U ^ oo and L — > oo (L the size of the system) do 
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not commute as is required for a hard projective model 
(no double occupancy in the original electron basis) to 
be the true low-energy theory of the Hubbard model. 

No such problem besets our low-energy theory. We can 
recover the original Hubbard model from our low-energy 
theory by simply integrating over ipi. Although this is 
not a sensible thing to do from a low energy perspective, 
it can be done exactly. To see how this happens, we 
rewrite Eq. (jlSp including the frequency dependence, 



ij \ {i,k) 

uj + U 



1 1] 



stpj - t hnj + (w + U)M 



- ^Tr In M 

P 

To simplify the expression, note that. 



(19) 



and 



~ 4i4T4,o-C",^Cj- |Cj- X = n,,_<,c|_^Cj>nj- _o(20) 

which comes from |A4yCj,|Cj^j^. The Lagrangian 

becomes 



L 




TrhiM 



(21) 



which yields the Hubbard model upon integration over 
Lpi (see the Appendix for details). 

We conclude from this analysis that our low-energy 
theory permits immediate correspondence with the orig- 
inal Hubbard model; that is, we have not lost any in- 
formation regarding the high energy scale, unlike projec- 
tive methods. All information regarding the high-energy 
scale is encoded into the emergent charge 2e bosonic ex- 
citation and its interactions. The IR physics will be de- 
termined by examining the low energy dynamics of the 
electrons/holes and (f. 



B. What this theory is not 

Rather than decoupling the on-site repulsion term, 
we derived our low-energy theory by exponentiating a 
(5— functional constraint on the heavy field, D. Nonethe- 
less, one might contemplate that our theory could be ob- 
tained by more traditional schemes, for example by some 
sort of Hubbard Stratonovich (HS) decoupling scheme. 
Since the interaction in the Hubbard model is entirely lo- 
cal, any decoupling by means of introducing an auxiliary 
field would yield only local interactions. The auxiliary 
field in Eq. (jlSp clearly generates non-local interac- 
tions as well as on-site interactions. Hence, Eq. 
cannot be obtained from a HS transformation. However, 
the non-local terms are scaled by t/U . Hence, it might 
still be maintained that the local terms dominate and 
I in fact that they could be obtained by some sort of HS 
transformation. Consider the identity, 



a e 



(22) 



which is true as long as A G M+ and A1A2 = —UX. The 
standard HS transformation assumes that A1A2 G K"'' and 
A2 = A^;. This necessarily leads only to the — C/ Hubbard 
model with X = C|Cj^. However, the constraint that the 
exponent on the right-hand side of Eq. be real can 
be relaxed^ in which case it can be applied to the +U 
Hubbard model as well. Nonetheless, this procedure will 
not yield the non-local terms in our low-energy theory. 
Further, it does not permit a clean identification of the 
field associated with the high-energy degrees of freedom. 
That is, there is no D-ficld in any version of Eq. (|22p . 
Consequently, this procedure is a non-starter for the con- 
struction of a proper low-energy theory. 



C. Electron doping 

For electron doping, the chemical potential jumps to 
the bottom of the UHB. Consequently, the degrees of 
freedom that lie far away from the chemical potential no 
longer correspond to double occupancy but rather double 
holes. To coarse grain these degrees of freedom, we ex- 
tend the Hilbert space in a similar way as the hole-doped 
theory, defining a new field D which will be constrained 
to describe the creation of double holes. Mathematically, 
all results in hole doping obtained in the previous section 
can be transformed to electron doping via a generalized 
particle- hole transformation (GPHT), namely, 



~* 



(23) 



where Q = (tt, tt) and Di is a fermion operator asso- 
ciated with double holes. The new bosonic field, (f, 
is the Lagrangian multiplier defined by the constraint, 
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ffcon = s9J2i 'fi{Di - 9cI t-cI^^) + h.c. According to Eq. 
([5|), an appropriate Lagrangian for the extended theory 
at electron doped can be constructed, 



L 



-t E 9i3 [A 
- {OCi^aVaCj-aDi + h.C.) 



(24) 



that preserves the distinct hops in the Hubbard model. 
Two differences to note are that 1) because the chemical 
potential resides in the UHB, the electron hopping term 
now involves sites that are at least singly occupied and 2) 
the order of the Di and Ci operators is important. If we 
integrate over (fii and then Di, all the unphysical states 
are removed and we obtain as before precisely Lnubb- 
Hence, both theories yield the Hubbard model in their 
UV limits. They differ, however, in the IR as can be seen 
by performing the integration over Di. The correspond- 
ing integral is again Gaussian and yields 



where 



u 



ij fj 



u 



sign change in the conserved charge will manifest itself 
as a sign change in the chemical potential as long as 
(iflifii) 7^ 0. Likewise, the correct — > oo limit is ob- 
tained as before. 



D. Half Filled Models 

At half filling when the chemical potential lies within 
the Mott gap, both double hole and double occupancy 
lie far away from the chemical potential. As two dif- 
ferent degree of freedom per site need to be coarse 
grained, we introduce two new fields D and D which 
when constrained correspond to double occupancy and 
double holes, respectively and extend the Hilbert space to 
(g)^ (JF| (X) JF^ (X) J^D ^ ^d)- The corresponding low-energy 
theory will be obtained by integrating over both D and 
D rather than by solving the constraint. As a result, in- 
tegrating over D and D do not yield identical results as 
would be the case if double occupancy or double holes 
were integrated over. 



1. Anderson Impurity 

To illustrate the process of coarse graining two high 
energy fields, we begin with a simpler model, the 
Anderson^ impurity model. 



(25) 



k,(T 



as the IR limit of the electron-doped theory. In Eq. ([25)) . 
the matrix Mij = Aiji is defined via the GPHT on the 
M matrix in Eq. pl|) . As i^i is a complex field, the GPHT 
interchanges the creation operators of opposite charge. 
We again make the identification s ^ t because the last 
term can also mediate spin exchange. 

When the boson vanishes, we do recover the exact 
particle-hole symmetric analogue of the hole-doped the- 
ory. Because the field now couples to double holes, 
the relevant creation operator has charge — 2e and the 
conserved charge is Q = J2i a 4 a''^i,o- - '^flfi- This 



where Oo- destroys an impurity electron and Ck.o- destroys 
a continuum electron. By setting Ef — —U/2, it costs an 
energy U/2 to create an double hole or a doubly occupied 
states on the impurity site which is analogous to the half- 
filled Hubbard model. In the following, we would like to 
show, by introducing two heavy fields D and D which 
correspond to the doubly occupied or double hole states 
on the impurity site respectively, the Kondo model with 
additional coupling to new bosonic fields can be derived 
perturbatively if Uk <C C/. 

The appropriate extended Hamiltonian is 



+sO(p'^ {D - Oai^ai) + h.c. + s9(p^ [d - da\a}^ + h.c. 



(26) 



In the current model, two bosonic field and (p are in- troduced which corresponding to the two constraints on 
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D and D fields respectively. If we first integrate out (/?, (p 
which result in delta functions corresponding to the re- 
moval of the unphysical states and then integrate out D 
and _D, Eq. (|25p is obtained precisely. This constitutes 
the ultra-violet (UV) limit of this theory. Similar to the 
previous discussions, we can derive the IR limit of the 
Anderson impurity model by first integrating out both 
D and D fields, which amounts to substituting 



D 



-0~|^V+E^kC^^-«^-j (27) 
I js>-K5]z;kC^y.an (28) 



into H . Wc finally obtain 
2 

HiR = i/Kondo - ■^(s^(^V + ■s^'^V) - (s'/'^ + #^)aTai 

- —{S(f^ + Sif>)'^V]^aaVaC\^^a + h.c. (29) 
k,cr 

where i?Kondo is the Kondo Hamiltonian 



U . 



ondo 



k,(T 



E 

k,k' 



4i;kfk' 
U 



(Ck'ffCk' ) • Si 



If. 
'2CkCk'(nT 



(30) 



Here, S^^p = a)^~^ a'^ is the spin operator of the impu- 
rity site. Thus the IR limit consists of the Kondo model 
in addition to the coupling between the electron and the 
new bosonic degree of freedom. 



2. Half-filled Hubbard Model 

Next, we perform the same procedure for the Hubbard 
model at half-filling by introducing two new fermionic 
fields D and D associated with the double occupancy and 
double holes, respectively. We consider the generalised 
Lagrangian 



L = 



t 

U 



')Ci„VcrCn-aDj + h.C. 



with the constraint terms given by 

Hcon = s9'^(pl{Di - 9ci^^Ci^i) + h.C. 

i 

+ ~s6Y,^\{b,-6cyi^) + h.c. (32) 

■i 

Here, (fi and Lpi are two the bosonic field with charge 
2e and — 2e respectively. Similar to the previous result, 
if we first integrate out both the bosonic fields Lpi, (pi 
and then Di , Di , the Hubbard model is obtained and the 
generalised theory Eq. [21] yields the correct UV limit. 
However, a different IR limit is obtained if we first inte- 
grate out Di and Di, 



H 



hf 
IR 



i i 

+ ^ ^(s<y9| + Sip,)h + h.C. 

i 

- ^ {s(pl - s(^i)ci,|Ci^| -I- h.C. 



E 2s^ 4. 2s^ ~t ~ \ 



(33) 



(31) 



This Hamiltonian is invariant under the transformation 
Ci,cr — > exp(iQ • Ri)c|j^, ipi ^ (pi and s <-> s. This invari- 
ant reflects the symmetry between the double occupancy 
and the double hole in the system at half filling. In con- 
trast to the doped case as in Eqs. and (|25p . no Ai 
matrices appear in the IR theory at half filling. Conse- 
quently, we arrive at a closed form for the low-energy the- 
ory at half-filling. The iD^b + bh^ terms include a spin-spin 
interaction as well as a three-site hopping term. However, 
at half-filling, the three-site hopping term vanishes. As a 
result, charge dynamics appear solely from motion of the 
charge 2e boson. This state of affairs obtains because at 
half-filling, charge dynamics persists only on the energy 
scale U. Since it is the boson that encodes the high- 
energy scale, it stands to reason that only the boson term 
mediates charge transport. As we show in Appendix A, 
it is the fi terms that break the spurious local SU(2) 
symmetry^ of the Heisenberg model and reinstate the 
global SU(2) symmetry of the Hubbard model. The lo- 
cal SU(2) symmetry in the Heisenberg model arose en- 
tirely from projection. This local SU(2) symmetry was 
noticed quite some time earlier-SS, but its origin was never 
clarified. In fact, it is straightforward to check that in 
the lower-Hubbard band at half-filling, all pcrturbative 
terms in t/U only mediate spin physics and hence pre- 
serve the local SU(2) symmetry of the Heisenberg model. 
Our derivation (see Appendix B) lays plain that as a re- 
sult of the boson terms, this symmetry is absent from the 
true low-energy theory of the Hubbard model. Note, in 
the Anderson impurity models, the bosonic terms play 
a similar role to destroy the local SU(2) symmetry that 
appears as a result of projection. 

Second, the true low-energy model preserves the sum 
rules associated with the original model. An essen- 
tial property^i^ of the half-filled Hubbard model when 
U ^ t is the presence of a surface in momentum space 
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(the Luttinger surface) where the single-particle Green 
function vanishes. Since the Mott state at half- filling has 
a gap, the non-trivial implication of the zero surface is 
that the real part of the Green function, 



i?.(k,0) 



de' 



°°d,'MSl£:)(34) 



vanishes. Here Ao-(fc,e) is the single-particle spectral 
function which we are assuming to have a gap of width 
2A symmetrically located about the chemical potential 
at e = 0. Because A(k, e) > away from the gap, and 
e changes sign above and below the gap, Eq. can 
pass through zero. For this state of affairs to obtain, the 
piecesof the integral below and above the gap must be 
retained. Projected models which throw away the UHB 
fail to recover the zero surface. What Eq. ([55)1 makes 
clear is that all the information regarding the surface of 
zeros is now encoded into the bosonic fields ipi and (pi . On 
the Luttinger surface, the self-energy diverges, represent- 
ing a break-down of perturbation theory. As the bosonic 
fields cannot be obtained from perturbation theory, we 
conclude that it is the emergence of the bosonic field that 
accounts for the breakdown of Luttinger's theorem'^^ and 
ultimately the Mott insulating state. 



E. Electron Operator 

In each of the low energy theories, the operator which 
creates a single electron represents a composite excita- 
tion. To determine its form, we add to each of the start- 
ing Lagrangians a source term that generates the canon- 
ical electron operator when the constraint is solved. The 
appropriate transformation that yields the canonical elec- 
tron operator in the UV, is 



L 



fl 



)4 



VrrD]eC, 



+ h.c. 



in the UHB for electron doping) with a renormaliza- 
tion from spin fluctuations (second term) and a new 
charge e excitation, Lp^^Mj^Ci-^. Consequently, we pre- 
dict that an electron at low energies is in a superposition 
of the standard LHB state (modified with spin fluctua- 
tions) and a new composite charge e state described by 
Ci-aM~^tf^y It is the presence of these two distinct ex- 
citations that preserves the dynamical (hopping depen- 
dent) part of the spectral weight transfer across the Mott 
gap. As shown in companion paper^^, there are also ex- 
perimental ramifications for the composite structure of 
the electron. 

At half-filling, a similar trick can be applied to generate 
the electron operator. In this case. 



L 



h.c. 



is the correct transformation to generate the canonical 
electron operator in the UV. If we now integrate the par- 
tition function over Dt and Di, we find that the electron 
creation operator at half-filling 



u 



has two important differences with its counterpart for 
n ^ 1. First, it lacks the standard LHB and UHB 
components as the chemical potential lies between both 
bands at half-filling. Second, the propagator Ai is ab- 
sent. Nonetheless, the electron at half-filling still has 
two components both above and below the chemical po- 
tential. The simplification that Ci^-afl (that is, the A4 
matrix is absent) constitutes the new charge e excitation 
may make subsequent calculations of the strength of the 
binding between the boson and a hole at least tractable 
within the framework of the Bethe-Saltpeter equations. 



III. FINAL REMARKS 



However, in the IR in which we only integrate over the 
heavy degree of freedom, 1?^, the electron creation oper- 
ator 



(35) 



For electron doping, we apply the generalized particle- 
hole transformation to obtain 



S 

U' 



t 



(36) 



as the generator of electron excitations in the IR. For ei- 
ther doping, the electron operator contains the standard 
term for motion in the LHB, (1 — ni _„)c\^ {ni-aCi^a) 



We have shown that a true low-energy theory of a 
doped Mott insulator possesses degrees of freedom which 
do not have the quantum numbers of the electron. The 
degree of freedom is a local non-retarded charge 2e bo- 
son and hence stands in stark contrast to the charge e 
boson in the slave bosoni2ii^i2i approach in which a di- 
rect integration of the high energy scale is not possible. 
Fundamental to theory here is that the boson does not 
act in its own Fock space, in contrast to other Fermi-Bose 
models'^^'^'^. That is, there are no free charge 2e boson 
states just as there are no free quark states in confin- 
ing theories. Rather, the charge 2e boson mediates new 
electronic states by forming composite excitations. As 
such the charge 2e degree of freedom is detectabl o^^i'^° 
through the substructure it provides in the electron ex- 
citation spectrum. In addition, the boson is not min- 
imally coupled to the electromagnetic gauge field. It 
acquires dynamics and hence a gauge coupling through 
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high order terms in the A4 matrix, essentially t^/U'^. At 
half-filling, the bosonic mode preserves the Luttinger sur- 
face on which the self-energy diverges or equivalently, the 
single-particle Green function vanishes. Since the boson 
represents a non-perturbative effect, it is not surprising 
that the Luttinger surface cannot occur without it. In a 
future publication^ we explore the role of the boson in 
mediating the normal state properties of the cuprates as 
well as the possibility that the Mott state is ultimately 
characterized by charge neutral bound states mediated 
by the hidden charge ±2e boson. 



review it to avoid any confusion in our derivation. To 
illustrate the method, we will begin with the familiar ex- 
ample of the non-linear cr-model, that is a bosonic theory 
with spherical target manifold. In Lorentzian signature, 
we introduce the spherical constraint by writing the cor- 
responding functional (5-function as an integral of a com- 
plex exponential, with Lagrange multiplier a: 



IV. APPENDIX A: LAGRANGE MULTIPLIERS 

Here we offer some details about the mechanics of La- 
grange multipliers. Though this is standard stuff, we 

I 



Zl[J] 



(37) 



r 



which after Wick rotation becomes 

Ze[J] = J[d(t)''da] e--/''^'^ 4(38) 

= J [da] e-^s=//('^) (39) 

with 

Seff (a) ^\Tr H-d^ + a)-^J^a (40) 



In eq. (|39p . we have performed the (jf functional in- 
tegrations. To proceed with the analysis of the model. 



J 



we investigate (j40l) by expanding a around its vev {a). 
It is crucial though to appreciate that if we go back to 
eq. ([38|) . we see that in order for cr to be a Lagrange 
multiplier field (in the Euclidean formulation), the fluc- 
tuations in cr should be taken along the imaginary axis 
in field space. That is, we write 



aix) = (a) H ^A(a;). 



(41) 



For uniform (cr), we then obtain 



Trln(-a2 + a) = -Tr In Ao + ^Ao(0) / X{x) + ^ f Ao(a;, y)A(y)Ao(y, x)A(x) + . . . 



(42) 



where Aq = —d" + (a). We thus find 



My,_iTrlnA 
2g 2 



r 



Ao(a;,y)A(y)Ao(y,x)A(a;) 



(43) 



Thus we see that there is a stable saddle point giving the 
familiar gap equation 



1 



d'^p 



(44) 



Now, in the theory considered in this paper, we have 
a similar situation. In the Lorentzian signature, we have 



Z ^ j[Vc Vc' VD VD'f V(p Vip^ exp^fd''x{Lo+H,,,^)^ 



(45) 
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where 



in 



d^0 



(1 



04 



3 



i,3,<y 



course in the reverse order, the Lpi integrals simply yield 
b— functions. In both cases, one ends up with the +J7 
Hubbard model. 



V. APPENDIX B: ABSENCE OF LOCAL SU(2) 
SYMMETRY 



(46) 



Here again, we have written the functional (5— function 
constraint as the integral of a complex exponential. In 
this case, however, ifi the Lagrange multiplier is a com- 
plex field such that (/?* = 3fi(/3i — V^ipi. After Wick rota- 
tion, we obtain the path integral in Euclidean signature. 



Z= [Vc Vc^ VD VD^ Vifi Vifi^] exp- ^''"="('^0+^-") 



(47) 

In either the Lorentz or Euclidean signatures, fi is a 
Lagrange multiplier. That is, the integral over ipi must 
still yield a functional 5— function even though the i is 
absent. This requirement dictates that we must integrate 
the fluctuations of both the real and imaginary parts of 
ipi along the imaginary axis as in Eq. (j4ip . The result 
is a stable Gaussian integral (if D is integrated first and 
then ifii) which can be evaluated using Eq. (I22p . Of 



Here we consider the possibility of trivially extending 
our model at half-filling to a locally SU(2) invariant the- 
ory. Let us organize the electron operators in the form 



(48) 



where the spatial index is suppressed. A local SU(2) 
transformation acts via left multiplication by an SU(2) 
matrix h so that 4* — > h'i>. Let us write h as 



h = 



a 



a* 



(49) 



where 



= 1. The electron bilinear C|C| is a mem- 



ber of the triplet of the local SU{2). We will now show 
that the correct "middle" term in the electron triplet is 
(n| + ni — 1)/a/2. Transforming the prospective triplet 
we find 



n-f +nj^ — 1 



V 



V2af3* 



V2a*f3 



*2 



a 



( CJ,C| 



V 



(50) 



By inspection it is clear that this matrix is unitary since h 
is unitary, and also that its determinant is (|Q;p-|-|/3p)^ = 
1. Thus, we have an SU(2) transformation and we con- 
clude that the electron triplet is 



W 



I 



\ 



n-j- 4-^1]. — 1 



(51) 



Thus, to make our theory local SU{2) invariant, we 
must also take D and ip to reside in triplets as well. Thus, 
we posit a charge zero bosonic field ipo which we collect 



into a triplet $ — {ip,(fo,(p)^ . SU{2) invariance would 
necessitate an additional constraint term of the form 



0{n^+ni-l)) (52) 



Just as in the cases of the earlier boson fields, the new 
field ifo enters the theory as a Lagrange multiplier cor- 
responding to a dynamical field Dq, which will be in- 
tegrated over. Now that we have determined the form 
of the electron triplet we can write the generalized La- 
grangian, now including the (po terms, in the form 
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PD]D, + -DlM 



9rj 



6c, 



+ 9DoiC, 



He, 



where the constraint is given by 

Hcon = se^ip\[D^ - ec^^^c^i) + h.c. + seY^(p\{D., - ecl-^c^^) + h.c. + sO'^ipliiDoi 



(53) 



(n,,| + - 1))(54) 



The undetermined coefRcients may be fixed by the condition that the theory reduces to the Hubbard model when the 
constraints are solved. Integrating over the tpi and then D fields, yields 



L 



E 



- b{l - ni^„)cl^Ci, 



(/? - + riii - 1) + {a - P + S)ni^n,i 



(55) 



up to total time derivatives. In order for this to yield the Hubbard model we must therefore have P — 1 — a and 
6 = 2(1 — a). We note that a = P = 25 = 1/2 solves these, and is in fact the SU{2) symmetric point. Similarly, we 
find a = —b = 1. The coefficient c is unconstrained because Dq has trivial dynamics. We have found then that the 
theory at half filling may be extended to an SU (2) invariant theory. This SU (2) acts only globally however; this may 
be plainly seen by examining the interaction terms on the second line of eq. I|53p . In order to make these terms local 
SU{2) invariant, an explicit SU{2) gauge field (a Wilson line) would have to be introduced. We conclude, then, that 
the local SU(2) symmetry of the Heisenberg model is broken by the presence of the bosonic degrees of freedom in our 
model. That a local SU{2) symmetric version of the theory cannot be constructed is not surprising as the Hubbard 
model lacks this symmetry - in that case it is broken by hopping terms (which again could only be made invariant 
by the introduction of explicit Wilson lines). In the case of the Heisenberg model, the local SU{2) symmetry appears 
strictly because of projection. Hence, the exact low-energy theory constructed without using projection should not 
possess symmetries not found in the Hubbard model. 
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